Intervortex quasiparticle tunneling and electronic structure of multi-vortex 
configurations in type-II superconductors. 



O 

o 

(N 
>^ 



(N 
(N 

o 
o 

u 

Oh' 



X3 
O 

o 



> 
oo 

in 
o 

o 



I 

o 
o 



X 
S3 



A. S. Mel'mkov, M. A. Silaev 
Institute for Physics of Microstructures, Russian Academy of Sciences, 603950, Nizhny Novgorod, GSP-105, Russia 

(Dated: February 6, 2008) 

The electronic spectrum of multi-vortex configurations in type-11 superconductors is studied taking 
account of the effect of quasiparticle tunneling between the vortex cores. The tunneling is responsible 
for the formation of strongly coupled quasiparticle states for intervortex distances a < a^, where the 
critical distance Oc is of the order of several coherence lengths 5. Analyzing the resulting spectra of 
vortex clusters bonded by the quasiparticle tunneling we find a transition from a set of degenerate 
Caroli - de Gennes - Matricon branches to anomalous branches similar to the ones in multi-quanta 
giant vortices. This spectrum transformation results in the oscillatory behavior of the density of 
states at the Fermi level as a function of a and could be observed in mesoscopic superconductors 
and disordered flux line arrays in the bulk systems. 

PACS numbers: 74.25.Jb, 74.25. Op, 74.78. Na 



Theoretical investigations of the quantum mechanics 
of low energy quasiparticle (QP) excitations in isolated 
vortices and vortex arrays are of crucial importance for 
understanding the basic thermodynamic and transport 
phenomena in the mixed state of type-II superconduc- 
tors at low temperatures. This issue has been studied 
for several decades starting from the pioneering work by 
Caroli, de Gennes and Matricon (CdGM) J^. For an iso- 
lated singly-quantized vortex there appears a so-called 
anomalous branch of the subgap spectrum which corre- 
sponds to the QP states bound to the vortex core because 
of the Andreev reflections from the gap potential profile 
A = |A(r)|e'^, where {r,9,z) is a cylindrical coordinate 
system, and z— axis is chosen parallel to the vortex line. 
The low energy spectrum of these localized states can be 
written as follows: E = — /iw ~ — /iAo/(fcj_^), where Ap 
is the gap value far from the vortex core, ^ is the su- 
perconducting coherence length, k± = ^Jk'^p — k1 is the 
Fermi momentum projection on the plane perpendicular 
to the vortex axis, and the angular momentum quantum 
number fx is half an odd integer. The overlapping of QP 
wave functions of neighboring vortices should perturb the 
CdGM spectrum and result, e.g., in the band structure 
effects in periodic vortex lattices 0- With the increase 
in the vortex concentration at high magnetic fields close 
to the upper critical field Hc2 there occurs a crossover to 
the Landau-type quantization for QPs precessing along 
the whole Fermi surface (see and references therein). 
A scenario of such crossover from discrete CdGM levels 
to the Landau-type spectrum is an appealing problem 
which would allow to understand the nature of de Haas - 
van Alphen oscillations observed experimentally even 
for magnetic fields H ~ (0.3 — 0A)Hc2 when the vortices 
are well-separated. The transformation of the QP spec- 
trum in the latter case should be controlled by the effect 
of QP tunneling between the vortex cores. This tunnel- 
ing phenomenon is expected to play an essential role also 
for exotic vortex configurations formed in mesoscopic su- 



perconductors of the size of several coherence lengths. 
In such systems the balance of competitive forces acting 
on vortices due to the screening current flowing at the 
sample boundary and the intervortex repulsion, leads to 
the formation of small size multi-vortex configurations 
(vortex molecules) and multi-quanta (giant) vortices 0|. 
These vortex states can transform into each other via 
magnetic field-driven first or second order phase tran- 
sitions. According to the general theory |^ the num- 
ber of anomalous branches (per spin) crossing the Fermi 
level at certain impact parameters fJ.F/k± ^ ^ equals to 
the winding number M of a multi-quanta vortex. The 
anomalous branch with = responsible for the peak 
in the density of states (DOS) at the vortex center exists 
only for vortices with an odd winding number. Gener- 
ally the spatial distribution of the DOS has the shape 
of rings with radii of the order of ^ 0i ■ The splitting 
of a multi-quanta vortex into M— vortex molecule with 
the intervortex distance a results in the transformation 
of these rings in the DOS profile into a set of peaks at 
the centers of individual vortices. The initial stage of 
such DOS transformation has been studied in Q by the 
perturbation method valid for a <JC ^. 

In this Letter we address the limit a > ^ of well- 
separated vortices in mesoscopic superconductors and 
disordered vortex arrays in the bulk systems well below 
Hc2- The goal of the present work is to study a scenario 
of the QP spectrum transformation caused by the for- 
mation of vortex clusters bonded by the QP tunneling 
between the Andreev wells in the cores. Note that here- 
after we neglect all the possible normal scattering effects 
at impurities or mesoscopic sample boundaries (lol 
which also can affect the low energy spectrum. 

Let us start with a qualitative analysis of the intervor- 
tex tunneling effect and consider a set of vortex lines 
parallel to the z— axis. In the plane (xy) the vor- 
tex centers defined as zeros of the superconducting or- 
der parameter are positioned at certain points r^. For 
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FIG. 1: Fig.l. Angular momentum /i as a function of 9p 
E = for (a) two and (b) three vortices 



at 



QPs propagating along the classical trajectories paral- 
lel to = k±{cos 9p, sin 9p) we introduce the angu- 



lar momenta /i — [r, k^] • zq = A;^rsin( 



6*) and 



fli = fj, — [vi, k^] • Zq defined with respect to the z— axis 
passing through the origin and with respect to the i-th 
vortex axis passing through the point , correspondingly. 
Neglecting the QP tunneling between the vortex cores we 
get degenerate CdGM energy branches: Ei = —Lojli. For 
a fixed energy E we can define a set of crossing branches 
on the plane (/i, fJ-i{Op) = —E/uj + [ri,kj^] • zq. 
These branches are shown in Fig ^ for two vortices with 
Ti = (—a/2,0) and r2 = (a/2,0) and three vortices at 
the apexes of the equilateral triangle with the center at 
the origin: ri = (0,a/A/3), ra = (a/2, -a/(2^/3)) and 
rs = (—a/2, —a/ (2-^3)). Each crossing point of branches 
fj-i{0p) and fJ.j{dp) corresponds to the trajectories pass- 
ing through the cores of i—th and j— th vortices. It is 
natural to expect that the degeneracy at these points 
will be removed if we take account of a finite probabil- 
ity of QP tunneling between the cores. The tunneling 
is determined by the exponentially small overlapping of 
wave functions localized near the cores and results in the 
splitting of energy levels: 6E ^ Aq exp{—kpaij /{k±£, j), 
where a^ = |rj — rj|. The estimate for the branch split- 
ting shown in Fig^Jby dash lines reads: 
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exp 



(1) 



As a result, we get the branches fi*(9p) with a quali- 
tatively new behavior: each of these branches consists 
of parts corresponding to the classical QP trajectories 
passing through the cores of different vortices. Keep- 
ing in mind general conditions of the semiclassical ap- 
proach validity it is natural to expect that a reasonable 



criterion on the intervortex tunneling efficiency can be 
obtained if we compare the splitting (5/x(aij) with the 
quantum mechanical uncertainty A/Lt of the angular mo- 
mentum. The latter value can be estimated from the 
uncertainty principle A^A9p ~ 1, where the expression 
for the angle uncertainty near the branch crossing points 
reads: A9p ~ Afi/{k±aij). Provided the branch split- 
ting is rather small, i.e. 5ii[aij) <C Aji ~ \/k±aij for 
all crossing points, the branches iJ.i{9p) are almost inde- 
pendent. In the opposite limit diJ,{aij) ^ \/k±aij the 
QP states of the z— th and j—th vortices appear to be 
strongly coupled by tunneling: near the branch crossing 
point the QP trajectory performing a precession in the 
core of the i—th vortex experiences a transition to the 
core of the j—th vortex. According to the above condi- 
tion on 5^{aij) the tunneling is most efficient for k± = kp 
and Qij < ttc, where ac — (C/2) hi[kpS^) is a critical inter- 
vortex distance. Using the percolation theory language 
we can consider the vortices to be bonded if a^ < ac 
and define a cluster in a disordered flux line system as 
a set of M vortices bonded either directly or via other 
vortices. The number M grows with the increase in the 
vortex concentration and becomes infinite at the percola- 
tion threshold or for a periodic vortex lattice. Certainly 
in mesoscopic superconductors the cluster dimensions L„ 
can not exceed the sample size R. Our further consider- 
ation of the bonded QP states is restricted to the limit of 
finite clusters with L„ <C tl, where ~ k±(j)Q/H is the 
cyclotron orbit radius, and (po = Trhc/\e\ is the fiux quan- 
tum. The cluster is characterized by a set of hybridized 
QP states; with a change in the k^ direction the wave 
function experiences a number of subsequent transitions 
between the cores of neighboring vortices. Taking, e.g., 
the upper branch in Fig^a) we obtain the wave function 
concentrated near the cores of the right and left vortices 
for the angular intervals < 9p < tt and n < 9p < 2tt, 
respectively. Further decrease in the intervortex distance 
results in the increase in the tunneling probability and, 
thus, the increase in Sfi{aij). Finally, for a.y —^ we get 
a set of M lines fj, — const parallel to the 9p axis, i.e. M 
anomalous branches crossing zero energy at angular in- 
dependent impact parameters and corresponding to the 
M— quanta vortex. Certainly this limit can be realized 
only in mesoscopic samples. 

Following ^2 one can obtain the discrete energy lev- 
els applying the Bohr-Sommerfeld quantization rule for 
canonically conjugate variables ^ and 9p: 



^l{9p)d9p = 27r(n-f /?), 



(2) 



where n and t < ne < M are integers, 27™^ is the period 
of the n{9p) function, and /? = 1/2 (/? = 0) for odd (even) 
ng values. Depending on the ratio Sfi{aij ) / A/x one should 
apply this quantization rule either to the branches ^ii{9p) 
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or to the branches fi*(6p). In the momentum region 



/cFy 1 - (min(aij)/ac)2 < Ifc^l < kp (3) 

one can neglect the branch sphtting S^{aij) ^ at all 
and Eq.lPl written for the branches iJ.i{9p) gives us the 
CdGM spectrum with a minigap cjq/S = ^^{kz = 0)/2. 
For min(aij) > Oc the CdGM expression holds for the 
entire momentum range. For vortices forming a cluster 
the QP states bonded by intervortex tunneling appear in 
a finite momentum interval 



\kz\<.kF\Jl- (min(ay)/ac)2 < kp ■ (4) 

In this limit the QP tunneling between the cores results 
in the qualitative modification of spectrum which can be 
obtained by substituting n*{Op) into Eq.Q: 

i?m(fc.)«y (^^ + 6.(ri,..rM)) , (5) 

where i = 1..M. The spectrum ISJ is similar to the one 
of a multi-quanta vortex 0, Qj ^] which recovers in the 
limit Uij when \bi\ < ^. The multi- vortex clus- 
ter geometry and its dimensions Ly determine the ef- 
fective impact parameters biiji, ..tm) which vary in the 
range —L^ < bi < L^. Taking a two- (three-) vor- 
tex molecule with ^ < a < Oc as an example we get 
bi.2 ~ ±a (6i,3 ~ ±a, 62 = 0). Contrary to the CdGM 
case the spectrum branches (jSJ can cross the Fermi level 
as we decrease the characteristic intervortex distance a 
and the minigap is suppressed. The DOS consists of 
M sets of van Hove singularities corresponding to the 
extrema of Eni{kz) branches. The energy interval be- 
tween the peaks belonging to each set is ujq. For a cer- 
tain fixed energy the DOS as a function of a exhibits 
oscillations with the period of the order of the atomic 
length scale Sa ~ kp^. Experimentally the intervortex 
distance can be controlled by a varying magnetic field. 
For typical values a ~ i/c^o/Jf we get the following field 
scale of DOS oscillations: SH/H ~ a/ Hloh /£f, where 
luh — \e\H/mc is the cyclotron frequency, m is the elec- 
tron effective mass, and ep is the Fermi energy. The 
oscillatory behavior should affect both thermodynamic 
and transport properties at low temperatures though 
in real experimental conditions the DOS peak structure 
is certainly smeared due to the various mechanisms of 
level broadening, e.g., finite temperature, fluctuations 
in vortex positions, impurity scattering effects, etc. It 
should be noted that for typical values kp^ — 10^ — 10'^ 
the critical distance Oc/C ~ 2 — 3 exceeds the core ra- 
dius and the spectrum transformation starts at the fields 
H ^ 4>o/(ic ~ Hc2[^n{kp£,)]^'^ when the vortices are in- 
deed well-separated. 

Now we proceed with a quantitative analysis of the 
intervortex tunneling on the basis of the Bogolubov - de 



Gennes theory: 

Hoi> + aa.ReAi' - o-j^ImA^ = , (6) 

where Ha = azUP - e&zA/cf ~ h^kl)/2m, * = iU,V), 
U and V are the particle- and hole - like parts of the 
QP wave function, ct, are Pauli matrices, p = — z/iV, 
A — A{x,y) is an order parameter profile. For extreme 
type-II superconductors we can assume the magnetic field 
H = —Hzq to be homogeneous on the spatial scale 
and take the gauge A = [H,r]/2. Within the quasi- 
classical approach the wave function in the momentum 
representation can be taken in the form: 

+00 

*(P) = -^ y" dse-*(IPl-'*'=^)^/'"V(s,0p) . (7) 

— 00 

The equation for ^(s, 9p) reads: H^p = Etp^ where 

H = -ihV^azg- + a^rReA- &ylmA+ ^^fi , (8) 

mV± — hk±, A = A{x,y), jl ~ —id/d6p, and the terms 
quadratic in H are neglected assuming Ly <C r^. The 
coordinate operator takes the form: 

r = .5ki/fci + {[ki,zo],/i}/(2fci) (9) 

where {A, B} = AB + BA is an anticommutator. Let 
us emphasize that the Hamiltonian (jSJ takes account of 
noncommutability of fl and 6p and, thus, the above de- 
scription involves the angular momentum quantization. 
Replacing /i by a classical variable we get Andreev equa- 
tions along straight trajectories. For an isolated vortex 
positioned at r = and described by the gap function 
Ai(r) = At,(r)e'^ a standard solution of the above equa- 
tions corresponding to low energy CdGM levels reads: 
Tp — ex.Y){i^9p)ijjQ{s,6p), where 



The physical picture of the intervortex QP tunneling 
effect can be illustrated by consideration of the simplest 
generic problem: two neighboring vortices positioned at 
Y± — (±a/2,0). The distances from this vortex pair to 
other vortices or to a mesoscopic sample boundary are 
assumed to be larger than the critical distance Cc so that 
their influence can be neglected. To describe the case of 
well-separated vortices a > ^ it is natural to use the tight 
binding approximation for the QP wave function: 

{l^ = c+{9p)f+Ms.0p) + c^{9p)f^Ms,0p) , (10) 
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where the operator 
1± = exp 



acos6'.„ , , 



ik I a 



transforms the CdGM states of a single vortex at the 
origin into the eigenstates for isolated vortices at r± = 
(±a/2, 0), ip± is the order parameter phase induced at the 
center of right (left) vortex due to the presence of neigh- 
boring vortices or boundaries. The coefficients c± satisfy 
the condition c±{9p + 2tt) = —c±{6p), so that the total 
wave function is single valued. In the angular intervals 
I sinflpi > £,/a the QP trajectories can not pass through 
both vortex cores and we have two solutions with either 
c+ = or c_ = 0, i.e. the states corresponding to one of 
the vortices. The only effect of neighboring vortices or 
boundaries on the QP energy for such trajectories could 
be associated with the Doppler shift caused by the su- 
perfluid velocity field induced by the external sources in 
a given vortex core. However this superflow would pro- 
duce the Lorentz force acting on a vortex and, thus, 
should vanish for a static vortex configuration. This con- 
dition gives us the relation a = a{H) which depends on 
the vortex arrangement and sample geometry (see, e.g.. 

The angular intervals | sin 6'p | < ^/a should be consid- 
ered separately since the corresponding QP trajectories 
pass through both vortex cores. Substituting Ea. Hl()(l 
into the quasiclassical equations with the Hamiltonian 

multiplying them by the functions ^T±^o^ i Sind in- 
tegrating over s, one obtains: 



i(7z fc^a(cos — 1) ; 



rC = Ec 



(11) 



where c = (c+,c_). For the trajectories with By close 
to or TT we find an approximate expression for the 
intervortex tunneling amplitude: t ~ Aqc""'^^/^^'^^^. 
Note that here we omit small Qp— dependent corrections 
to w resulting from the angular dependent distortions 
of the gap profiles in the cores and neglect the renor- 
malization of the w and t values caused by the terms 
cx hjjii LOolS^/a)'^ <C coq. Let us introduce the function 
b — exp{iE9p/uj — ik±aaz{cos6p — l)/2)c. For \6p\ ^ ^/a 
the Ea. lfTT)) reads: 



-ujfib - 



Opa^b 



tdxb — 



(12) 



This problem is equivalent to the one describing the inter- 
band tunneling [l3j or the one-dimensional (ID) motion 
of a Dirac particle in a uniform electric field and the so- 
lution can be written in terms of the parabolic cylinder 
functions £',vv(cr): 



6 = 



(ilAa(CT) + d2D^ai-a) 



y/-ia{diDia-i{(T) - d2Dia-i{-cr)) 



(13) 



where a — t"^ /{uj^k±a) ~ (J^(a)/A^)^, and a — 
Opy^k±a/i. Taking the asymptotical expressions for 




FIG. 2: Fig.2. Two- vortex molecule {kp^ = 200): (a) en- 
ergy spectrum for a = 2.1^; (b) the DOS per spin projection 
averaged over the interval 5E = O.OSojq around i? = as a 
function of the intervortex distance a. The CdGM spectrum 
is shown by the dash lines, uo = fcF/(4ajo) is the averaged 
CdGM DOS 



Dia{<j) in the interval max[Q;, l\/^/k±a <C \0p\ < ^/a we 
get: b{0p > 0) = Sb{9p < 0), where 



((TyRer + (Tajlmr) e 



2F(1 - ia) _^a/2 • w \ 
r = T{a) = — e ''"''^ smh(7ra) , 

\'—2TTia 

I is the unity matrix, 7 = k± a9^ /A + ahi{\dp\\/k±a), and 
r(a::) is the gamma function. For the angles 9p close to 
TT the solution can be obtained using a transformation: 
b{\ep - 7r| < (/a) = axb{\ep\ < ^/a). Matching the 
wave function in various angular domains and using the 
condition c{6p) — ^c{0p + 2tt) we obtain the spectrum: 

cos(7r£:/w) = ±e~''"/2 ^3 sinh(7ra) sin(7rx) , (14) 

where nx ~ k±a + a\n{k±£^^ /a) + arg(r(l — ia)) + tt/A. 
The generalization of such matching procedure for M— 
vortex clusters is straightforward. The spectrum calcu- 
lated using the Eg. 1)14(1 is shown in Fig[21[a) for a typical 
parameter kp^ — 200 and for a rather small intervortex 
distance a < Qc- One can see that the transformation of 
the spectrum E{kz) occurs according to the scenario sug- 
gested above: as we decrease the distance a below Oc the 
crossover to the double quanta vortex spectrum starts 
in the region of small fc^ values defined by the condition 
a > 1. For a <C 1 we get the CdGM spectrum with a 
small oscillatory correction: 



E-uj{n+l/2) ~ =F(-l)"t^V2a/7rsin(fc±a + 7r/4) . (15) 



The effective minigap Emin = wo(l/2— \f2afK) vanishes 
for a ~ Cc. For intermediate values \Jk\^^ ja ^ a ^ 1 
the spectrum takes the form Q if we put 61.2 = ±x/fc^. 
This expression can be matched with the one obtained 
using a standard WKB expansion for Eq. Hll|) in the limit 
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a 1 if we put 



In agreement with general arguments presented above 
the DOS in the hmit a < Gc consists of two sets of peaks 
shifted by the value ujq{2x — [2x]), where [...] denotes the 
integer part. The oscillatory behavior of the DOS at the 
Fermi level as a function of a is shown in Figl^Jb). 

In conclusion we note that the bending of vortex lines 
in a cluster can strongly affect the behavior of the wave 
functions along the magnetic field direction. Let us take, 
e.g., a two- vortex molecule and assume the function a(z) 
to change slowly on a scale ~ fc^^. One can replace the 
momentum in Eq.(|SJ) by the operator = —id/dz 
and consider the spectrum branches as a set of effec- 
tive Hamiltonians Eni{kz) describing the ID quantum 
mechanics of QPs. Such approach is analogous to the 
one suggested in |l4j for QPs in superconducting/normal 
metal layered structures. Taking a simple model profile 
a{z) = a — az"^ / and choosing n ~ ^kpa/ir such ID 
problem for energies close to the Eni{kz) extremum can 
be reduced to the one of a quantum mechanical harmonic 
oscillator with the effective mass to* ~ rnkpS,^ /a and in- 
terlevel spacing fl ^ uJoV^A/ L. The corresponding wave 
functions are localized at a length scale ~ L{kpdS^)~^^^ . 
Such localization effect should result, in particular, in the 
suppression of the heat transport in the z— direction. 
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